Stability of Bose-Einstein condensates in a circular array 
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The properties of the superfluid phase of ultra cold bosonic atoms loaded in a circular array 
are investigated in the framework of the Bose-Hubbard model and the Bogolhibov theory. We 
derive and solve the Gross-Pitaevskii equation of the model to find that the atoms condense in 
states of well-defined quasimomentum. A detailed analysis of the coupling structure in the effective 
quadratic grand-canonical Hamiltonian shows that only pairs of distinct and identical quasimomenta 
are coupled. Solving the corresponding Bogoliubov-de Gennes equations we see that each pair 
of distinct quasimomenta gives raise to doublets in the excitation energy spectrum and that the 
quasimomenta of the zero-energy mode and of the occupied state in the condensates are identical. 
The dynamical and energetic stabilities of the condensates are determined by studying the behavior 
of the elementary excitations in the control parameters space. Our investigation establishes that 
superflow condensates exists only in the central region of the first Brillouin zone whereas there is none 
in the last quarters since they are energetically unstable, independently of the control parameters. 

PACS numbers: 03.75.Fi, 05.30.Jp 



I. INTRODUCTION 

The experimental realization of Bose-Einstein conden- 
sates in a periodic potential created by interference of 
laser beams, called an optical lattice [USUI, opened up 
the possibility of studying the properties of ultra cold 
atoms in a variety of conditions. Examples are the ac- 
cess to strongly interacting atomic regimes [3, @ and to 
states which are fundamental for applications in quan- 
tum information Q. These systems have proven to be a 
rich research field that offer the unique possibility to in- 
vestigate fundamental questions from condensed matter 
physics to quantum optics. 

A Bose-Einstein condensate in an optical lattice is 
nearly a perfect physical realization of the Bose-Hubbard 
model, a fact first pointed out by Jaksch and coworkers 
■ This model gives a description of the physics of inter- 
acting bosonic atoms trapped in a lattice potential with 
atomic tunneling between nearest neighbors and on-site 
repulsion. The main advantage of these systems is that 
they are highly controllable. As an example, when the 
ratio between the interaction and hopping terms is varied 
by changing the laser intensity, a quantum phase transi- 
tion from a superfluid phase to a Mott insulator phase is 
observed ||. 

In this paper we use the Bose-Hubbard model and 
the Bogoliubov theory to study the dynamical and en- 
ergetic stabilities of Bose-Einstein condensates loaded in 
a periodic ring. The dynamical and energetic stability 
of condensates in optical lattices have been investigated 
both theoretically EH and experimentally 

[II EE EE EE El ■ The purpose of this paper is to give 
a complete and self-contained description of the prop- 
erties of condensates, its excitation spectrum and stabil- 
ity. Briefly, we characterize the condensates by the quasi- 
momenta of occupied states which reveals the existence 
of equilibrium current carrying condensates [IE EH , we 
analyze the coupling structure of the effective quadratic 



grand-canonical Hamiltonian in the quasimomenta basis 
which allowed us to establish the composition of the ele- 
mentary excitations, the doublet structure of excitation 
energy branch and the proper assignment of the zero- 
energy mode [l9| and the phonon limit. We find that in 
the zero-energy mode, which is a direct consequence of 
the atom number violation by the Bogoliubov theory, the 
quasimomenta of the atoms are identical to the quasimo- 
mentum of the occupied state in the condensate. Turning 
to the phonon limit we find that it is achieved when the 
relative quasimomentum I — p goes to zero, with I and 
p being, respectively, the quasimomenta of the excita- 
tion and of the occupied state in the condensate, a fact 
overlooked in the literature [IE HE] ■ 

The stability of the condensates is determined by the 
behavior of the excitation energies and composition of the 
elementary excitations in the control parameters space. 
In fact an equilibrium state is dynamically stable if all 
excitation energies are real, [IE EH E2I; 13 1 , and it is 

nil ei mi. 
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energetically stable if they are positive, 
In our study we identify two mechanisms of energetic 
instability: "crossing" that occurs when the excitation 
energy vanishes and change its sign and "no-crossing" 
that occurs when the excitation energy is strictly nega- 
tive, independently of the control parameters. From this 
analysis we determine the domains in the control param- 
eters space where it is possible to find metastability in 
the system, that is, dynamically and energetically stable 
condensates. Metastable current carrying condensates 
correspond to local minima of the energy and they are 
candidates to present superfluid motion [111 [l8| . 

This paper is organized as follows. In the Section II 
we derive and solve the Gross-Pitaevskii equation for the 
Bose-Hubbard model to show that the atoms condense 
in states with well-defined quasimomentum [Io| . reveal- 
ing the existence of equilibrium current carrying states 
[18j |. Besides, these states with well-defined quasimomen- 
tum define a single-particle basis which diagonalizes the 
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hopping term of the Bose-Hubbard Hamiltonian. In the 
Section III we determine the energies and the composi- 
tion of the elementary excitations by the diagonalization 
of the effective grand-canonical Hamiltonian of the Bo- 
goliubov theory. We express it in the quasimomcntum 
representation to show that only pairs of quasimomcnta 
are coupled. We identify these pairs to cast the effective 
grand-canonical Hamiltonian into the form of a sum of 
terms each one involving pairs of identical and distinct 
quasimomenta. This fact allow us to reduce the process 
of diagonalization of a 2M x 2M matrix, with M be- 
ing the number of lattice sites, to one of 2 x 2 matrices 
when the quasimomenta of the pairs are identical and 
4x4 matrices when they are distinct, which implies that 
the effective grand-canonical Hamiltonian is diagonal in 
blocks. In fact, when M is odd, we have M ~ 1 pairs with 
distinct quasimomenta and only one pair with identical. 
On the other hand, for M even, we have 4p — 1 pairs with 
distinct quasimomenta and two pairs with identical. We 
will see that 2x2 blocks correspond to one excitation en- 
ergy whereas 4x4 blocks correspond to doublets. In the 
Section IV we determine the dynamical and energetic sta- 
bility of the condensates by studying the behavior of ex- 
citation spectrum and the composition of the elementary 
excitations in the control parameters space. We found 
that the stability properties of the condensates depend 
only on a combination of control parameters r = ^j, 
where n = j^, J and U are, respectively, the density of 
atoms, the hopping and the on-site strenghts. We deter- 
mine critical values of r that define, in the control pa- 
rameters space, the domains where it is possible to find 
superflow states in the system. The dynamical and en- 
ergetic phase diagrams are shown. A summary and our 
conclusions are presented in the Section V. 



the number of atoms on the Ath lattice site. These opera- 
tors satisfy the periodic boundary condition a\ = ax+M- 
The first term in Bose-Hubbard Hamiltonian is the hop- 
ping term that describes the tunnelling of atoms among 
neighboring lattice sites with hopping strength J and its 
effect is to delocalize the atoms over the lattice. The 
second term describes the inter-atomic on-site repulsion 
with interaction strength U whose effect is to localize the 
atoms on the sites. 



B. The condensates 

In a condensate all atoms are in the same single- 
particle state and therefore the many-body state can be 
written as 



pfiV 



(2) 



where ft is the creation operator of the state occupied 
by the atoms 



M-l 



ft 



-y 



z x a\ 



A=0 



which requires that the complex parameters z\ must sat- 
isfy the constraint 



M-l 

E 

A=0 



N. 



(3) 



II. THE GROSS-PITAEVSKII EQUATION OF 
THE BOSE-HUBBARD MODEL 



The Bose-Hubbard Hamiltonian 



The parameters z\ are determined by a variational 
principle in which we minimize the mean value of Hbh 
in the boson condensate \c), given by ((2]), subject to the 
constraint j3| . The minimization lead us to the equation 



The physics of an ultra cold, dilute and interacting 
Bose gas in a latticepotential is captured by the Bose- 
Hubbard model @, HiJ. This model can be seen as an 
one-mode approximation that involves the states in the 
first band of the optical lattice. In the tight-binding ap- 
proximation, the homogeneous Bose-Hubbard Hamilto- 
nian for a system of bosons in a periodic circular array 
with M sites is given by [1(3] 



- J(z A +i + za-i) + U - —j \z x \ 2 z x = (J,z x (4) 

which is the Gross-Pitaevskii equation for the Bose- 
Hubbard model with /i being the chemical potential. A 
solution of this equation that satisfy the constraint (J3]) 
has the general form 



H 



BH 



M-l 

-'£( 

A=0 



u 



M-l 



a{ax + i+a{ +1 ax) + - ^ 

A=0 



h\{h\ - 1). 

(1) 

In (fTJ), a\ and at are, respectively, bosonic annihilation 
and creation operators of atoms on the Ath lattice site 
and fi\ = alax is the atom number operator that counts 



and, from ((4]), the chemical potential takes the explicit 
form 



M= -2Jcos0+^(l-i 



(5) 
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The possible values of 9 are fixed by the periodic bound- 
ary condition z\ + m = z\ which requires that (e* 6 ') = 
1. This equation has M solutions which are the Mth 
roots of unity, that is, 9 k = "^ff where k is an inte- 
ger defined on the set { — 4f + 1, ■ ■ ■ , 4f} for M even or 

(M-l) (M- 
2 ' " " " ' 2 



for M odd. It follows from all 



these considerations that the bosonic atoms condense in 
states with well-defined quasimomentum 



4- 



>M 



M-l 
A=0 



These states form a single-particle basis which diagonal- 
izes the hopping term of the Bose-Hubbard Hamiltonian. 
Thus, in this representation, the Hamiltonian (TT|) has the 
form 



Hbh = £ e kA{A k 



+ Mfci +k 2 -k 3 - /c 4 )4 4 i t fe3 i fe2 i fcl (7) 



2M 



{fc 3 } 



where 



e k = — 2 J cos 



2nk 
W 



and the Kronecker modular delta 8m (k) is equal to one 
if k is an integer multiple of M and zero otherwise. 



III. THE ELEMENTARY EXCITATIONS 

A. The effective grand-canonical Hamiltonian and 
the coupling structure of quasimomenta 

To derive the Bogoliubov-de Gennes equations we first 
define the shifted operators c k and cl by 



Ai 



Zq&k,q ~t~ ^k : 



where z q is the condensate wave function with quasimo- 
mentum ^jj- . Next we write the zero temperature grand- 
canonical Hamiltonian f2 = Hbh — I^N , with N denoting 
the number operator, as a normal order expansion with 
respect to the shifted operators, c&, 



(8) 



state of the operator When we calculate the mean 
value of Q in this vacuum state it is clear that the only 
contribution comes from the term luq whose minimization 
with respect to z q lead us to equation 



-2 J cos 



2nq U . 

— 1 H 

M M 1 



(9) 



Besides if we calculate the mean value of N in this ground 



state, we get N = 



Since this equation plus 



determine z q except by a phase factor, it can be taken 
(6) as a real parameter. From these considerations it follows 



that z n 



N and 



2irq N TT 
/, = -2Jcos— + -U. 



(10) 



Notice that the expressions (JS|) and PH| , for the chemical 
potential, coincide in the limit N ^> 1. 

In the Bogoliubov theory the main hypothesis is that 
the majority of atoms remains in the boson conden- 
sate state and thus the dynamics of the system in the 
neighborhood of equilibrium states is described by the 
quadratic term <1>2 since the equation ([9]) makes Cj\ iden- 
tically zero. Thus, neglecting the third and fourth or- 
der terms in (J5J) , we obtain the effective grand-canonical 
Hamiltonian 



nU 



6 M(k + k'- 2<?)(44 + c k ,c k ) (11) 



k,k' 



where 



M) 



efc 



2J cos 



27rq 



2irk\ 
M J 



(12) 



The diagonalization of the effective Hamiltonian deter- 
mines the energies and the composition of the elementary 
excitations. The non-diagonal term in (|11[) is the one 
that creates and annihilates pairs of atoms whose quasi- 
momenta are fixed by the condition <5a/ (fc + k' — 2q) = 1 
which is satisfied if (fc + k' — 2q) is an integer multiple of 
M, 



k + k' - 2q = v (q) M. 



(13a) 



As the quasimomenta are restricted to the first Brillouin 
zone, the possible values of v (q) are given by 



where the normal ordered operators d), involve i shifted 
operators. 

In the Bogoliubov theory the ground state is given by 
the vacuum of the shifted operators which is a coherent 



for M odd 

v{q) = I 1 0, 1 for M even 
.0, -sgng 



if q = 
if q^O 



(13b) 
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where sgaq is the signal function which is equal to 1 
if q > and -1 if q < 0. In the u(q) = equation 
both quasimomenta are in the first Brillouin zone whereas 
when v{q) ^ one of the quasimomenta is outside, the 
effect of v(q) ^ is to bring it back to the first Brillouin 
zone. Notice that, except in the case q = and M odd, 
the pairs of quasimomenta coupled in (TTTJ) are determined 
by two equations. 

From these equations two properties of the coupling 
structure of quasimomenta emerge. The first one is the 
absence of ramifications, that is, one quasimomentum 
appears only once in the set of pairs satisfying (|13a|l . 
The second one is that the sets of pairs of quasimomenta 
obtained from these two equations are disjoint. These 
considerations indicate that when M is odd, we have 
M ~ 1 pairs with distinct quasimomenta and only one pair, 
(q, q), with identical. On the other hand, for M even, we 
have M _ i 

pairs with distinct quasimomenta and two 
pairs with identical, (q,q) and (q + ^§-v(q),q + ^-f(q)j, 
with v{q) = — 1 if q > and v(q) = 1 if q < 0. In 
the Tables U and HH we identify the pairs of quasimo- 
menta coupled in pip for M even and odd, respectively, 
irrespective if they differ by an exchange of the quasimo- 
menta. 



B. The Bogoliubov-de Gennes equations and 
elementary excitations 

The identification of the pairs of quasimomenta cou- 
pled in the effective Hamiltonian shows that it is block 
diagonal: 2x2 blocks when the quasimomenta of the 
pairs are identical and 4x4 blocks when they are dis- 
tinct. In fact, (jlip can be cast into the form 

(k,k) (*,*') 

k^k' 

where the first sum involves pairs with identical quasi- 
momenta while the second distinct, with fi2(k,k) and 
1i2(k,k') given explicitly by 



the main properties of the eigenvalue problem. In what 
follows we will discuss separately the diagonalization of 
blocks that involve identical and distinct quasimomenta. 



• Identical quasimomenta: 2x2 blocks 

When M is odd there is only one pair with identical 
quasimomenta, (q, q), which is equal to the quasimomen- 
tum of the occupied state. Solving the corresponding 
Bogoliubov-de Gennes equations (|Alj) for h,2(q, q), where 
A = B = nU, we find that the eigenmode has zero en- 
ergy and zero norm. This result is well known and it is 
a consequence of the broken continuous symmetry which 
arises due to violation of the atom number conservation 
by the Bogoliubov theory. In this case fi2(q,q) can be 
written as 



h2{q, q) 



EL 

2x 



nU 



where the hermitian operator p is given by + c q ) 

and the inertial parameter \ is equal to (2nU) . 

On the other hand, when M is even, we found 
two pairs with identical quasimomenta (q, q) and 

(q+fv(q),q+fv(<l)), wi th K?) = -1 if « > and 
v(q) = 1 if q < 0. The pair with quasimomentum equal 
to the occupied state, (q, g), gives raise to a zero-energy 
eigenmode as stated before. In the case of the other 
pair, (g + 4£i/(<7), q + ^v(q)) , the Bogoliubov-de Gennes 
equations are given by 



t t 27TQ 

4 J cos -jf 
nU 



nU 



nU 
4 J cos ^ 



nU 



9+^(9) 



1 

-1 



V q+M-v(q) 



(14) 



where we have used e 



(9) 



= 4Jcos^, from (TT 



and (|13ap . Solving these equations we find the excitation 
energy 



fi2{k, k) 



and 



(4 9) 



h 2 (k,k>) - +n U)c\c k + {e^ +nU)c\,c k , 
+ nU(c\c\, + c k 'C k ). 

The diagonalization of quadratic Hamiltonians of bosonic 
operators is a standard problem fully explained in refer- 
ence [l9| . In the Appendix [X] we discuss briefly the pro- 
cedure of this reference to fix the notation and to present 



nUy k c k + —{c' k cl 



(9) 



CfeCfe) 



E (i) 



1 4 J cos 



2tt<7 

IF 



4 J cos 



2irq 



2nU 



and the eigenvectors associated to 



E 



-E 



(8) 



V 



(9) 



and 



the 

(9) 



(15) 



pair 



respectively, with V^^ v{q) 



q+^ V (qy 9+^K«)V 9+f-^(9) ' 9+#f(g) ' 

( u i+f Hi) v i+f'Ai) ) i 
where the symbol T denotes the transpose of the 
corresponding matrix. As will be clear latter on, another 
useful quantity is the ratio of the components of the 
eigenvector given by 
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Table I: The set of pairs of quasimomenta coupled in the effective grand-canonical Hamiltonian for M even. 



M even 


q < 


/ jn ^ — 

u(q)=0: (f -2M-l,-f + l) , . 

i/(«) = l: -2|g|) (- 


.,(-M,-M),...,(-f + l,f -2|g|-l) 
g | + f,-kl + f),...,(f -2 k |,f) 


, = 


/ M 1 A* _i_ 1 \ 
1/(0) = 0: l^" 1 '-— + 1J 

1/(0) = 1 : 


,...,(0,0),..., 

{ 2 > 2 J 


(-f + l,f -1) 


q > 


/ Af M , „\ 
u(q)=0: {—'-—+ 2 1) 

u(q) = -l: ^_M +2? _ 1) _M +1 ^ ) ... j 


,...,(<?,,),..., 
(g-f.^-f) 


(-f + 2,,f) 
,...,(-# + l,-f +2,-1) 



Table II: The set of pairs of quasimomenta coupled in the effective grand-canonical Hamiltonian for M odd. 



M odd 


q < 


( M- 

u(q) = : V~ 
u{q) = 1 : 


±-2|,|,-^l) ,...,(-1,1,-1 
(^i,^fi-2|,| + l),..., 


7 |),...,(-A^i,A^- 
^-2| 9 | + 1,^1) 


%l) 


, = 


i/(0) = : 


(M_i,_M_i),... i( o,0),..., 


' M-1 M-l\ 

v 2 ' 2 ; 


q>0 


i/( 9 ) = : 
!/(,) = -1 : ( 


(^,-« f l+2,),...,(,,,) 
'-^1+2,-1,-^1),..., 




;-^+2,,^i) 


) 



u 



4 J cos ^ 



nU 



4 J cos 2g2 (4 J cos ^? + 2nf7) 



(16) 



• Distinct quasimomenta: 4x4 blocks 

In this case the Bogoliubov-de Gennes equations for 
h 2 (k, k') is given by 



/ e{ q) + nU 


nU 



-k 









+ n<7 


nU 


n?7 


4 9) + 









/ 1 
10 
0-10 

\ -1 



Tit/ 






nU j 





(17) 



It is easily seen that equation (fTT|) leaves uncoupled 
the vectors v} q) = (u k v k ,) T and V$ = 

(0 uy Vk ) T . Therefore this equation can be reduced 
into two 2x2 eigenvalues equations 



life/ 



e^ y ' +nU J \ v k 



= E\ 



(9) ( 1 o 



" V -1 



Wfc' 



(18b) 



Solving these equations we find that the excitation ener- 
gies are given by 



p(q)± _ Jg) 



(k, k') ± yje { f (k,k') (e^ (k, k>) + 2?i?7 



(19a) 



p(?)± _ J?) 



where 



(*, k') ± JJ§ ] (jfc, jfc') (e^ (fc, jfc') + 2nU 

(19 b) 



(/) 



nU 



nU nU 
e («) 

A:' 



(</) 



Ufc' 

1 

-1 



Vk' 



(18a) 



.(9) _ ,(?) 
e fc' 



£ a"* (^j k ) — 



with ej: 9) given by |T2 



e (9) + 6 (9) 



e^(fc,fc') 



Concerning the amplitudes, we found that they sat- 



isfy the following relations u 



(±) _ Jt) „,(±) _ „,(+) 



J k< > 
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MlT 1 ^ = and ul^ — vW* . These properties allow 

us to identify two excitation eigenmodes whose pairs 
of opposite energies are ^£^ 9)+ , — E^ + = E^f^ and 

(^Ejf + , —Ejf + = E < j q ' > ~^j . The corresponding eigenvectors 

are V^ q)+ , 7 ^ )+ = v}f~ and V% )+ , 7 V# )+ = , 

whose components are Vj ? ' + = ( v& ) and 



Thus, the diagonalization of 4 x 4 blocks gives raise 
to doublets of excitation energies which are degenerate 
when q — and q = 4p, the last one only for M even. 

One consequence of the above considerations is that 
the eigenmodes belonging to a doublet have the same 
norm, the ratio of the amplitudes given by 



,(+) 



,(+) 



.(9) 



(k,k')+nU 



^4 9) (M')(* 



(9) 



(k,k') + 2nU 



nU 



(20) 



IV. STABILITY ANALYSIS 



A. Dynamical stability 

According to the Bogoliubov theory an equilibrium 
state is dynamically stable if all the excitation energies 
are real. The existence of at least one complex energy 
is sufficient to guarantee the dynamical instability of the 
corresponding condensate. In what follows we investi- 
gate the dynamical stability of the condensates. We will 
analyze separately the cases where the quasimomenta of 
the pair are identical and when they are distinct. 

• Pair of identical quasimomenta 

As seen in the previous section, for both M odd and even, 
there is a zero-energy eigenmode with quasimomentum 
equal to the quasimomentum of the occupied state in 
the condensate. This is a consequence of the violation 
of atoms number conservation introduced by the Bogoli- 
ubov approach. This mode is always present and leads 
to an indifferent equilibrium which does not affect the 
stability of condensates. 

For M even there is one more pair of identical quasi- 
momenta equal to k — q + 4f z/(g), where u{q) = —1 if 
q > and v{q) = 1 if q < 0. From (JTSJ) it follows that 
the excitation energy is real if one of the conditions 



2?rg „ 

cos > 

M 



or 



27rg nU 



M 



2 J 



< 



(21a) 



(21b) 



is obeyed. The condition (|21a[) is satisfied by the conden- 
sates whose quasimomentum of the occupied state is in 
the interval < \p\ < |, with p = and the condition 
(|21bp is satisfied when the quasimomentum is in interval 
§ < \p\ < 1" such that 



r < 



2nq 



(22) 



where r is the combination of the control parameters . 

Up to now our analysis have established the domain of 
stability of a particular eigenmode k = q + ^-v{q). To 
establish the dynamical stability of the condensate we 
have to analyze the behavior of all the doublets in the 
control parameters space. 

• Pair of distinct quasimomenta: the doublets 



From (|19a|) and (|19b[) , the excitation energies of the dou- 
blets corresponding to the pair [k, k') are real if one of 
the two conditions 



e ( s q \k,k') >0 

or 

e (q \k,k') + 2nU <0 
is satisfied, where Eg (k, k') is given explicitly by 



(k,k'(k)) = 4Jcos 



2nq 



sin 



7r(/c — q) 
M 



(23a) 
(23b) 

(24) 



with k — q assuming the values shown in Table El 

From ([24| it is easily seen that the condition (|23a|) 
is satisfied for condensates with quasimomenta defined 
in the interval < \p\ < |. These condensates whose 
quasimomentum of the occupied state is in the central 
region of the first Brillouin zone are always dynamically 
stable, independently of the values of the control param- 
eters. On the other hand the condition (|23b[) is satisfied 
for condensates with § < \p\ < n that obey the inequality 



2irq . 2 ir(k - q) 
r < — cos — — sin — — . 

M M 



(25) 



7 



This condition refers to stability of a particular doublet, 
(k,k'(k)). However we are interested in establishing a 
condition that guarantee that the excitation energies of 
all doublets are real. This requirement is fulfilled if we 
take the minimum of the right hand side of inequality (|25|) 
which, according to Table W\ is achieved when k—q = ±1. 
Thus, it follows that the condensates whose quasimomen- 
tum of the occupied state is in the last quarters of the 
first Brillouin zone, § < |p| < 7T, are dynamically stable 
if they obey the condition 



2irq 2 7r 
r < — cos — — sin — . 

M M 



(26) 



Notice that the stability of the doublets implies the sta- 
bility of the equal quasimomentum pair since the inequal- 
ity (gll) is contained in ([55]). 

In the Table IIIII we summarize our findings and a dy- 
namical stability phase diagram is shown in Figure [1] for 
M = 5 sites. We disregard the condensates with quasi- 
momentum \p\ = ^ j the reason being that we cannot find 
a Bogoliubov transformation that diagonalizes the effec- 
tive grand-canonical Hamiltonian since all eigenvectors 
of the Bogoliubov-de Gennes equation have zero norm. 



Table III: Dynamical stability of the condensates. 



Dynamical Stability 


< |p| < f 


always stable 


| < |p| < 7T 


stable if r < — cos p sin^ 



condensates with the smallest \p\ in the set | < |p| < ir 
start to become unstable until r > sin 2 jt, when all the 
condensates in the last quarters of the first Brillouin zone 
are unstable. 

Recent experiments [lj| have shown that the onset of 
dynamical instability occur when two atoms in the con- 
densate can elastically scatter into a final state where 
they have quasimomenta different from p = -^9-. To 
point out the relevance of our analysis to this matter, 
we cast the stability condition (|23ap into the form lll| 



2e„ < e k + e k > 



(27) 



This condition reveals that in the condensates whose 
quasimomentum of the occupied state is in the central 
region of the first Brillouin zone, pairs of atoms can- 
not elastically scatter into two different quasimomentum 



states 



2irk 



'- and 



2irk' 



M M ■ 

For the condensates whose quasimomentum of the oc- 
cupied state is in the last quarters of the first Brillouin 
zone the stability condition (|23bj) takes the form 



2e„ > e k + ey + 2nU. 



(28) 



Notice that in this case the quantity 2e q — e& — e^/ is 
positive and has a lower limit equal to 2nU which implies 
that two atoms can inelastically scatter with the excess 
of energy being transferred to the other atoms in the 
system. However, notice that in the thermodynamical 
limit all these condensates are dynamically unstable. 



B. Energetic stability 




-0.8 -0.4 0.4 0.8 



P_ 
— 

Figure 1: Dynamical stability phase diagram for M = 5 sites. 
Notice that p assumes only discrete values 0, ±0.47r, ±0.87r 
indicated in the figure. The condensates with quasimomenta 
and ±0.47r, which are in the central region of the first Brillouin 
zone, are always dynamically stable. On the other hand the 
dynamical stability of the condensates with quasimomenta 
±0.87r, which are in the last quarters of the first Brillouin 
zone, depends on the control parameters. The continuous 
curve is only to guide the eye. 

When r = all the condensates are dynamically sta- 
ble. Actually they are eigenstates of the Bose-Hubbard 
Hamiltonian when U = 0. When r starts to increase the 



In the framework of the Bogoliubov theory an equilib- 
rium state is energetically unstable if there is at least one 
negative excitation energy. Recall that to an elementary 
excitation we associate a pair of eigenvectors with oppo- 
site eigenvalues (E, —E) and opposite sign of the norm, 
where the value of the excitation energy is the eigenvalue 
of the eigenvector with a positive norm. 

We can distinguish two mechanisms of energetic in- 
stability. One is "crossing" that occurs when a positive 
excitation energy vanishes and changes the sign and the 
second is "no-crossing" where we always have a nega- 
tive excitation energy. The difference between these two 
mechanisms is that "crossing" depends on the control 
parameters whereas "no-crossing" does not. 

We restrict our analysis of energetic stability to dy- 
namically stable condensates whose excitation energies 
are all real, a necessary condition of energetic stability. 
Our procedure to establish the energetic stability of the 
condensates is, for each eigenmode, to identify the eigen- 
vector with positive norm and the sign of the correspond- 
ing eigenvalue. The sign of the norm depends on the size 
of the ratio ~ . If < | ~ | < 1 the norm is positive 
and if | ^ | > 1 the norm is negative. Following what we 
have done in the case of dynamical stability, we will dis- 
cuss separately the case of identical and distinct pairs of 
quasimomenta. 



8 



• Pair of identical quasimomenta 

As pointed out before the zero-energy eigenmode does 
not affect the stability of the condensates. Therefore we 
are left with the pair k — q + 4pi/(q), with v(q) = — 1 if 
q > and v(q) = 1 if q < 0, that exists only for M even. 
By inspection of the equation (fT6|) we find that 



However we are interested in a condition that guaran- 
tee the stability of all the doublets. This requirement is 
satisfied if we take the maximum of the right hand side, 
which is achieved at k — q = ±1. Thus, the condensates 
such that < |p| < S are energetically stable if they 
obey the condition 



a) if cos ^ > 0, then V 



(?) 



has a positive norm; 



b) if cos ^ 



^ < 0, then V 



(?) 

?+¥"(<?) 



has negative 



2nq 



2irq 



cos " —rr + T cos — — > cos — , 
M M M 



Since in both cases the corresponding energy S^m,,^ 

is always positive we conclude that in case a) , where < 
|p| < -|, this particular eigenmode is always energetically 
stable. On the other hand, in case b), where § < |p| < n, 
it is always energetically unstable. 

• Pair of distinct quasimomenta: the doublets 

Inspection of equation (|20p shows that, 



which can be written in a more convenient way as 



/ 2nq 7r \ / 2irq n \ 2irq 

sm — — — sin — — - + — < r cos — — . 

\M M J \M M J M 



(30) 



a) if 4 9) 0, k 1 ) > 0, then and Vfi n have positive 

norm; 



rig)- 



b) if 4 9) (fc, k') + 2nU < 0, then V^ q>+ and V^ >+ have 
negative norm. 

Notice from (|19al) and (|19b[) that the sign of the en- 



(9) + 



r(«)+ 



ergies Ej 9 ^ and Ej 9 '^ depends on the term e^'ik^k'), 
given explicitly by 



(«)/ 



e%\kMk)) = 2J S ^s* 2 ^pl (29) 

for all the pairs of quasimomenta (fc, k'(k)). According to 
Table [V] we see that this term is semi-negative definite, 



Thus, in this case, the route to energetic instability is the 
"crossing" mechanism. 

In the case b) where the condensates are such that 
f < |p| < 7r, one of the energies of the doublets are al- 
ways negative. To see this, note that Ejj is always 

positive and the corresponding eigenvector, , has 

negative norm. Thus, by the Bogoliubov criteria, the ex- 



citation energy of this mode is —E 



(<?)+ = E (i)~ 



ii 



which 



is strictly negative, independently of the control parame- 
ters. Therefore, we conclude that the condensates whose 
quasimomentum of the occupied state is in the last quar- 
ters of the first Brillouin zone are always energetically 
unstable, in this case through the "no-crossing" mecha- 
nism. In the Table HVl we summarize our findings. 



e { 2 } {k,k'{k)) <0. 



This property implies that e\ Q j + is always positive, 

whereas Ej q ^ + can change its sign. 

In the case a), where p is defined in the range < 
\p\ < -|, the condition for the (k, k'(k)) mode to be ener- 
getically stable is E^ + > which can be cast into the 
form 



4 9) (k, k') I (k, k') + 2nU\- e ( %> (k,k r ) 



.(<?) 



.(<?), 



> 0. 



Inserting the expressions (f24l) and (|29f in the above in- 
equality we obtain 



2 2nq 2nq 2 ir(k - q) 
cos — — + r cos — — > cos — — . 

M M M 



Table IV: Energetic stability of condensates. 



Energetic Stability 


o < |p| < f 


stable if sin (p — -fj) sin (p + jj) < r cosp 


f < |p| < 7T 


always unstable 



In the Figure [2] we present the energetic stability phase 
diagram for M = 5 sites and in the thermodynamical 
limit, M 3> 1, n = -jg fixed. We plot the curve that de- 
fines the boundary of energetic stability only for r > 0. 
This curve is confined in the central region of the first 
Brillouin zone where the condensates are always dynam- 
ically stable, independently of the control parameters. 
In the last quarters, the condensates are always energet- 
ically unstable. 
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(a) 




-0.5 0.5 

P 



Figure 2: Energetic stability phase diagram (a) for M — 5 
sites and (b) in the thermodynamical limit M S> 1, n = & 
fixed. In (a), since r is a semi-positive definite parameter, we 
did not plot the curve that defines the boundary of energetic 
stability in the interval < \p\ < where r is negative. 
Notice that, in (a), p assumes only discrete values 0, ±0.47r, 
±0.8-7r indicated in the figure and the continuous curve is only 
to guide the eye. The two vertical continuous lines define the 
region where the condensates are always dynamically stable, 
independently of control parameters. 



One question that we can address is to determine 
the interval of control parameters in which there are 
metastable current carrying condensates, that is, dynam- 
ically and energetically stable condensates where the oc- 
cupied states have a finite quasimomentum. These su- 
perflow states correspond to local minima of the energy 
and they are candidates to present superfluid motion. 
To find this interval notice that we can define a critical 
value r, say r st , such that for r > r st the condensates 
with < \p\ < | are all metastable. If p max is the high- 
est value of p defined in the interval < \p\ < j, from 
(|30p the critical value r s t is given explicitly by 



= sin ( Pmax - §) sin (p max + ^) ^ 

COS Pmax 

We can also define another critical value r, say r unst , such 
that for r < r unst there is no metastability. When this 
occur only the p = condensate is stable. From this 
consideration follows that 



_ sin (p min - 1^) sin (p min + Jf) 

^unst — • 
COS p min 

Since |p min | = |f , r unst is given by 

sin (U) sin (A) , 
cos§- 

When r increases continuously from zero there is a hier- 
archy in the appearance of these superflow states, start- 
ing at r = r unst where pairs of degenerate condensates 
with quasimomentum p and —p, beginning with p = ^ , 
become metastable up to r > r st when all states are 
metastable. 



V. SUMMARY AND CONCLUSIONS 

In this paper we use the Bose-Hubbard model and the 
Bogoliubov theory to investigate the properties of ultra 
cold bosonic atoms loaded in a periodic ring with M sites. 
First we derive and solve the Gross-Pitaevskii equation 
of the model and from the analysis of the solutions we 
show that the atoms condense in states with well-defined 
quasimomentum whose values are the Mth roots of unit, 
restricted to the first Brillouin zone. Thus, besides the 
usual zero quasimomentum condensate, we have equilib- 
rium states with non-zero quasimomentum which corre- 
spond to current carrying condensates [181 ] . These states 
with well-defined quasimomentum form a basis that di- 
agonalize the hopping term of the Bose-Hubbard Hamil- 
tonian. 

Following the Bogoliubov theory we derive the effective 
grand-canonical Hamiltonian, quadratic in the shifted 
operators, whose diagonalization gives the energies and 
the composition of the elementary excitations. A detailed 
analysis of the coupling structure in the effective grand- 
canonical Hamiltonian shows that only pairs of identical 
and distinct quasimomenta are coupled. An immediate 
consequence of this coupling structure is that the effec- 
tive Hamiltonian is block diagonal: 2x2 blocks when the 
quasimomenta of the pairs are identical and 4x4 blocks 
when they are distinct. The diagonalization of the 2x2 
blocks shows that the pair (g, q) with quasimomentum 
equal to the quasimomentum of the occupied state in 
the condensate gives raise to the zero-energy eigenmode. 
On the other, the diagonalization of the 4x4 blocks 
gives raise to doublets of excitation energies which are 
degenerate when q = and, for M even, q = 4p This 
shows that when q ^ the excitation spectrum has a two- 
branch structure. We have also found, by inspection of 
the excitation energies (|19a[) and (|19b[) . that the phonon 
limit is achieved when the relative quasimomentum I — p 
goes to zero, with I — ^f^- being the quasimomentum 
of the excitation. Indeed, in this limit it follows that 
Eiji — cui\l — p\ with two different sound velocities: 
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c\ = 2 J(— sinp+^/r cosp) andc^} = 2J(sinp+^/r cosp). 
These properties are signatures of the finite size of the 
quasimomentum of the occupied state in the condensates. 

Our stability analysis shows that the condensates in 
the central region of the first Brillouin zone, < \p\ < 
^ , are always dynamically stable whereas the dynamical 
stability of the condensates in the last quarters, f < |p| < 
7r, depends on the control parameters. When r increases 
from zero, these condensates start to become unstable 
beginning with the one of smallest \p\ ending up with the 
instability of all condensates, when r > sin 2 jfe. 

Concerning the energetic stability, we found that the 
condensates in the last quarters of the first Brillouin zone 
are always unstable whereas the energetic stability of the 
condensates in the central region depends on the control 
parameters. We show that when r < r un st only the p = 
condensate is stable. However, when r > r unst , pairs of 
degenerate condensates with non-zero quasimomentum, 
p and — p, start to become energetically stable ending up 
with all the condensates in the central region of the first 
Brillouin zone being energetically stable, when r > r st . 

As discussed in the paper, metastable current carrying 
condensates are candidates to present superfluid motion. 
A metastable state is both dynamically and energetically 
stable, consequently a local minimum of the energy. Our 
analysis shows that the Bogoliubov theory predicts that 
there is an interval in the control parameters space where 
metastable current carrying condensates exist. The num- 
ber of these states increases with r and they are restricted 
to the central region of the first Brillouin zone. 
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Appendix A: THE DIAGONALIZATION OF A 
QUADRATIC HAMILTONIAN OF BOSONIC 
OPERATORS 

A general quadratic Hamiltonian of bosonic operators 
is given by 

i,j id 

where cj and Ci are creation and annihilation bosonic 
operators, Aij and Bij are elements of an hermitian and 
symmetric matrices, respectively. The diagonalization of 
this quadratic Hamiltonian consists in finding a canonical 
transformation to quasiparticle operators, 

&n = £«4-«?Cr). 

^ T 

such that H, when written in terms of this operators, 
takes the form of a system of non-interacting quasiparti- 
cles. From this constraint it follows that 



which is the matrix form of Bogoliubov-de Gennes equa- 
tions that determines the excitation energies E n and 
the composition of the elementary excitations V n — 
( u n v n ) , with T denoting the transpose of the cor- 
responding matrix. When the eigenvalues E n are real, 
they appear in pairs (E n , —E n ) with opposite sign of the 
norm, (\ u r\ 2 ~ \ v r\ 2 )- I n fact, if V n is an eigenvector 
with the eigenvalue E n then 7V™*, with 




is an eigenvector with the opposite eigenvalue, — E n , and 
opposite sign of the norm. The excitation energy is iden- 
tified with the eigenvalue whose eigenvector has a positive 
norm. 



Appendix B: PARAMETRIZATION OF THE 
PAIRS OF DISTINCT QUASIMOMENTA 

As seen before the pairs of distinct quasimomenta, 
(fe, fc'), coupled in the effective grand-canonical Hamil- 
tonian identify the doublets that compose the excita- 
tion spectrum of the condensates. However k and k' 
are not independent since they are related by equation 
(|13ap . Therefore we need one parameter to identify the 
doublets. Our choice was the relative quasimomentum, 
— q), which leads to an ordered and one-to-one 
parametrization of the doublets. This can be easily seen 
noticing that we can cast the equation Q13ap into the form 



(k-q) + (k' -q)=v (q) M, 

and, from it, follows the parametrization shown in the 
Table EO 

Table V: Identification of the doublets in the excitation spec- 
trum by the relative quasimomenta (k — q) of the atoms in 
the condensate. 



Condensate 


Doublets (k — q) 


AI odd 


M even 


q < 


Kk-q< 


l<fc-?<4f-l 


q > 


<k-q< -1 


-4f + l<fc-g<-l 



This parametrization can be extended to include the 
pairs of identical quasimomenta. Indeed, k — q = 
for the {q,q) pair and k — q = v{q)^ for the 
(? + !/(?)# ,q + v(q)f) pair. 
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